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Detection of binary FM by multiplication of the received signal by itself 
delayed is analyzed. Error rates vs signal-to-noise ratio for additive Gaus- 
sian noise are calculated as a function of sampling time, differential delay 
at the receiver, and delay distortion in the channel. It is found that the 
differential detector can give better performance than the more conventional 
zero-crossing counter or frequency discriminator under conditions of severe 
delay distortion in the channel. 

I. INTRODUCTION 

It has been found possible to realize excellent practical performance 
in FM transmission of binary data by use of a detector in which the 
signal is multiplied by a delayed replica of itself. This method has been 
called "differential detection" on account of its resemblance to the 
scheme of that name in widespread use as a detector of phase-modulated 
waves. The name "product demodulation" has also been applied. We 
can regard the detector either as a particular kind of frequency dis- 
criminator or as a phase comparator operating on the phase changes 
inherent in an FM wave. The former concept is suitable for a steady- 
state analysis, while the latter is more convenient in the study of signal 
transitions. 

II. THE DIFFERENTIAL DELAY PRODUCTOR AS A FREQUENCY DISCRIMINATOR 

Viewed as a discriminator, the detector has a steady-state response 
function calculable by multiplying a sine wave A cos (co c + w)t by the 
corresponding delayed wave A cos [(a> c + u)(t — r)], i.e., 

A cos (co e + tojt'A cos [(oi c -j- w)(t — t)} 

= (A 2 /2) cos (w fl + u )t (1) 

+ (A 2 /2) cos [2(« e + u)t - (a> c + «)t]. 
Ill 
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When the double-frequency component is suppressed by a low-pass 
filter, we obtain the response 

Vi S = (A 2 /2) cos (a, + «)r. (2) 

Consider w as the frequeDcy deviation in radians/sec from the midband 
frequency co c . Then a discriminator characteristic can be realized by 
setting cos co c t = 0, sin co c t = ±1, giving 

V lf = T(A s /2) sin cot. (3) 

The resulting steady-state response is nearly proportional to frequency 
deviation over the range in which sin cot is approximately equal to cot. 

A linear relationship is not necessary for binarj' FM detection, since 
only the sign of the deviation is significant. Unambiguous results in 
the noise-free case can be secured over a range of cot from -ir to it. 
In particular, if co = ±co 6 /2 where co 6 is 2t times the bit rate, values of 
t in the range zero to the reciprocal of the bit rate could be used, and a 
value equal to half the bit interval appears to be a good compromise. 
As t is made small the linearity improves, and, as will be shown later, 
the performance approaches that of an ideal phase differentiator, which 
we shall refer to as a <fo/dt detector. The latter type is of particular 
interest because its performance is closely approximated by either a 
zero-crossing counter or a frequency discriminator. 

III. FM OETECTION BY DIFFERENTIAL PHASE COMPARISON 

We illustrate the operation as a differential phase comparator by 
following a particular noise-free sequence through the detection process. 
The binary signal to be transmitted is shown in Fig. 1 (a) . It is assumed 
that this rectangular wave modulates the frequency of a carrier with a 
total shift between mark and space equal to the bit rate. This results 
in the phase change during a marking bit interval differing from that 
during a spacing bit interval by 360°. With respect to the mid-frequency 
as a reference, the variation of carrier phase versus time becomes ±180° 
per bit interval, as indicated by the solid triangular wave of Fig. 1(b). 
When the channel is shaped to give a raised-cosine pulse spectrum at 
the demodulator input, the phase-versus-time pattern becomes rounded 
at the transitions approximately as shown by the dotted waveform of 
Fig. 1(b). The received signal is passed through a network with an 
envelope delay of one-half bit interval and with a phase shift at midband 
frequency of 270°. The phase-versus-time pattern of the delayed signal 
is shown by the dashed-line wave of Fig. 1(b). For a long mark interval 
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Fig. 1 — (a) Binary data, (b) phase of undelayed and delayed FM waves, 
(c) demodulated signal, (d) signal space diagram. 
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the two phase patterns are seen to be in phase, while for a long space 
interval they are 180° apart. The demodulation process consists of 
taking the product of these two waves with a switch-type modulator 
and filtering out the high-frequency component. If the amplitude 
modulation of the received signal is negligible, the demodulated signal 
becomes equal to the cosine of the phase difference between the delayed 
and undelayed versions of the signal. 

The phase difference between the two signals at numerous points is 
indicated in Fig. 1(b). The shape of the demodulated wave, together 
with appropriate sampling times to recover the binary information, is 
shown in Fig. 1(c). It will be noted that the demodulated wave is de- 
layed one-quarter bit interval from the instantaneous frequency of the 
received signal but is advanced one-quarter bit interval from the in- 
stantaneous frequency of the delayed signal. 

A signal space diagram, such as described by J. R. Davey, 1 is given in 
Fig. 1(d) with lettered points corresponding to the lettered points on 
the phase pattern of the received signal. At a transition in the demodu- 
lated signal, the vectors representing the received signal and its delayed 
version are in phase at points such as (i), (k), (n), and (o). At a sampling 
instant having a transition on only one side, such as at t\ , the vectors 
representing the received signal and the half -bit delayed signal are at 
points such as (d) and (i). These form an angle of approximately 80° 
which, when the delayed signal is shifted 90°, becomes an angle of 10° 
or 170° depending on whether a mark or space is received. At a sampling 
instant having a transition on both sides, such as at t 2 , the two vectors 
are at points such as (o) and (n). The angle is then about 70°, which, 
after the added 90° shift, results in an angle of 20° or 160°. 

For a constant-amplitude signal the departures from the ideal 0° 
and 180° angles would decrease the amplitude of the demodulated signal 
at the sampling time. The diagram of Fig. 1 (d) shows, however, that at 
least one of the vector amplitudes at these times is above steady state. 
If the increased signal amplitude appears at the linear input of the 
demodulating productor, it more than compensates for the phase error, 
thus tending toward an overshoot of the baseband signal. The low -pass 
filtering will of course determine the final extent, if any, of the overshoot. 
It is of interest to note that a linear productor would over-emphasize 
these amplitude variations. Consequently, it would be expected that a 
switched-type modulator would result in a more perfect eye pattern.* 

* By "eye pattern" is meant the oscilloscope trace obtained by sweeping the 
detector output against a linear time base synchronized with the bit rate. A basic 
description of the properties of such patterns has been given by Brand and Car- 
ter. 2 
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This does not mean, however, that the probability of decision error 
would be different. In fact, as will be shown later, limiting one of the 
two inputs to the multiplier, and hence obtaining in effect a switched- 
type modulator, does not change the error rate. 

Fig. 2 shows a computer printout of the noise-free eye pattern corre- 
sponding to one-half bit delay when the demodulator follows a product 
law, the total frequency shift is equal to the bit rate, and full raised- 
cosine spectra apply. The origin is taken at the midpoint of the bit 
interval in the undelayed wave. Traces of like polarity are concurrent at 
this point, showing the absence of intersymbol interference at these 
particular instants. The peak responses of the detector are not reached 
until a time later by half the differential delay, and the individual peaks 




TIME IN BIT INTERVALS 



Fig. 2 — Eye pattern from product of undelayed and delayed FM waves with 
delay of one-half bit interval. 

for different adjacent data sequences are spread over a range of values. 
The traces concur again at a time equal to the differential delay, and 
the same signal levels are observed as at the first point of concurrence. 
The decision-threshold (zero signal level) crossings are almost concurrent 
at a time preceding the origin by half the difference between the bit 
interval and the differential delay. 

The time relations of Fig. 2 are peculiar to the choice of differential 
delay, ratio of frequency shift to bit rate, and spectral shaping. It will 
be shown later* for example, that if the differential delay is 5 bit in- 
tervals, the first concurrence of traces is at the origin, the second one is 
at 8 bit intervals later, the peaks occur at 5/2, and the threshold cross- 
ings nearly at —(1 — 8)/2. If the raised-cosine shaping were changed to 
any other satisfying Nyquist's first criterion 5 for suppression of inter- 

* See Figs. 9-12 and the discussion following equation (9'J). 
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symbol interference, all concurrences would be maintained except the 
threshold crossings, which could spread out more. Finally, if the rela- 
tion between frequency shift and bit rate were changed, all concurrences 
would be destroyed, but this would not by itself imply a significant 
increase in error rate. 

Eye patterns obtained when a limiter is inserted in one, but not both, 
of the two inputs are shown in Fig. 3. In Fig. 3(a), the undelayed signal 
is limited, giving in effect a switched modulator with the polarity of the 
undelayed signal switching the delayed signal. In Fig. 3(b), the delayed 
signal is limited, thereby interchanging the switching roles of Fig. 3(a). 
Comparing Fig. 3(a) with Fig. 2, we note that the concurrence of traces 
at 0.5 is destroyed and that the peaks of the responses are shifted to the 
left. In Fig. 3(b) the concurrence at is destroyed and the response 
peaks are shifted to the right. In spite of these differences, which can 
be verified by fairly straightforward analysis, the probability of error 
at a specified sampling instant must be the same for all three cases, as 
will be shown in detail later. 
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Fig. 3 — Eye patterns when one signal is limited before application to mul- 
tiplier input. Delay of one-half bit interval, (a) Undelayed signal is limited; (b) 
delayed signal is limited. 
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The eye pattern for a dtp/dt receiver, which is equivalent to the case 
of zero differential delay, is shown in Fig. 4. The traces concur at the 
peak response, which is at the origin, and the threshold crossings are 
half a bit interval away, as would be deduced by setting 8 = in the 
discussion of the detector with differential delay. As 5 approaches zero, 
we thus approach an equivalence with the performance of conventional 
frequency detectors. 

The question of utility of eye patterns for nonlinear detection proc- 
esses merits some further discussion. Although we cannot deduce error 
rates from them, we can at the least distinguish between "go" and "no 
go" conditions in the absence of noise. We can also use a given pattern 
as a basis for choosing the best sampling time. In making the choice, 
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Fig. 4 — Eye pattern for dip/dl detector. 

we consider both the horizontal and vertical margins, since relative 
immunity to timing jitter is important as well as the spread of sample 
values at the time of decision. We have verified in our calculations that 
the direction of change in the error rate with sampling time can be 
deduced by comparing eye openings for the two instants. Such deduc- 
tions are not valid for comparing error rates corresponding to eye 
patterns for two different conditions. 

Another use of the eye pattern is in laboratory diagnosis of system 
distortion. For such purposes we make use of established correspondence 
between the nature of the eye and various kinds of distortion in the 
particular system under test. 



IV. THE EVALUATION PROBLEM 

We are concerned here with the performance of the differential FM 
detector when the binary FM signal suffers transmission impairments. 
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We analyze particularly the effects of delay distortion and noise in the 
channel. A significant measure of performance is the curve of error rate 
vs signal-to-noise ratio at a specified bit rate for a channel with specified 
amplitude and delay variations with frequency. It is important that 
such a measure be determined over the range of transmission impair- 
ments encountered in actual channels. It is found that the principal 
virtue of the differential-delay scheme relative to the more familiar 
axis-crossing and frequency-discriminator types is an improved im- 
munity to severe delay distortion. 

Since there are many parameters which influence the performance of a 
data transmission system, the discussion would get out of hand if 
individual attention were given to all possible combinations of condi- 
tions. Fortunately, we can select representative regions of interest 
which are describable in terms of relatively few quantities. Our approach 
makes use of both direct analysis and digital computer back-up. In 
calculating the noise-free responses of the system under various condi- 
tions, we first establish formulas for a generalized data sequence. We 
then go to the computer to evaluate response functions vs time for a 
family of sequences of given length. 

In theory it would be possible to calculate error rates by adding 
programmed noise in the computer calculation of response functions. 
A practical deterrent to such a procedure is the long computing time 
needed in the interesting cases of almost error-free transmission. It has 
been our experience that determination of error rates by computer 
simulation of additive noise is inferior to computer evaluation of analytic 
formulas for probability of error. Use of the latter procedure has been 
chiefly successful with additive Gaussian noise, but this does not neces- 
sarily imply serious limitation of utility. The premise that rank-order 
established on the basis of Gaussian noise holds for other kinds of 
interference has a good empirical foundation. A useful simplification 
from the Gaussian analysis is that the curves for error rate vs signal-to- 
noise ratio tend to be roughly parallel for different systems and to be 
characterized sufficiently well by their asymptotic slopes. 

In terms of the normalized signal-to-noise ratio M, which is defined 
as the ratio of average signal power to the average white Gaussian noise 
power in a bandwidth equal to the bit rate, the asymptotic error rate is 
expressible 3 in the form F(M) exp (—kM). The function F(M) turns out 
to be of slight interest because of its relatively minor effect when M is 
large. For practical purposes the number k determines the performance. 
For an ideal binary phase modulation system, k has its maximum value 
of unity. The quantity 10 logio (1/k) expresses noise impairment in db 
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relative to the ideal. It has been found possible to calculate this quantity 
directly without determining the entire error-rate curve. 

The calculations we have made evaluate the effect of the following 
factors important in system design: 

1 . The sampling time relative to the signaling interval. The preferred 
sampling time is indicated by the eye patterns, but is more precisely 
established by error-rate calculations. 

2. The length of the delay line. Equivalence of zero delay with a 
flip/dt detector establishes a reference at one end of the range in terms of 
a better-known system. 

3. The data sequence. Results for the most and least vulnerable 
sequences are exhibited. 

4. The delay distortion. Parabolic and linear variation with frequency 
are studied. The results are presented in terms of maximum delay 
variation expressed in bit intervals. 

V. THE MODEL 

A block diagram of the transmission system under study is shown in 
Fig. 5. The data source emits a sequence of binary symbols which for 
full information rate are independent of each other and have equal 
probability. The analysis can be generalized without analytical incon- 
venience by assigning a probability mi to one of the two binary symbols 
and 1 — mi to the other. In conventional binary notation the symbols 
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Fig. 5 — Binary FM transmission system with differential delay detection. 
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are 1 and 0. It is convenient to express binary frequency modulation of 
an oscillator in terms of positive and negative frequency deviations. 
The combination of data source and low-pass filter is accordingly defined 
by the shaped baseband data wave train 

•(*) = E b n g(t - nT) (4) 

71= — 00 

where 

b n = 2a n - 1. (5) 

The values of a„ represent the data sequence in binary notation. The 
probability is Wi that the typical a n is unity, and 1 — Wi that it is 
zero. The value of b„ is +1 if a n is unity, and —1 if a n is zero. The 
function g(t) represents a standard pulse emitted by the low-pass 
filter for a signal element centered at t = 0. 

Ideally, the oscillator frequency follows the baseband signal wave 
s(t). This would imply an output voltage from the FM oscillator speci- 
fied by 



V(t) = A cos \u c t + d + m / s(\)d\ . 



(6) 



Here, A is the carrier amplitude, oj c is the frequency of the oscillator 
with no modulating signal applied, / is an arbitrary reference time, 
0o is the phase &tt = to, and n is a conversion factor relating frequency 
displacement to baseband signal voltage. The instantaneous angular 
frequency of the wave (6) is defined as the derivative of the argument 
of the cosine function. It can be written in the form co c + w,- , where 
a),- , the deviation from midband, is ideally expressed by 

w» = iis(t)- (7) 

In the practical case, the transmitting bandpass filter restricts the 
frequency-modulated wave to the range of frequencies passed by the 
channel. The purpose of this filter is to prevent both waste of trans- 
mitted power in components which will not reach the receiver and 
contamination of the line at frequencies assigned to other channels. 
The result is a transformation of the voltage wave (6) to a band- 
limited form, which must depart in more or less degree from the ideal 
conditions of constant amplitude and of linear relationship between 
frequency and baseband signal. The line also inserts variations in ampli- 
tude- and phase-versus-frequency which cause further departures from 
the ideal. For our purposes it is sufficient to combine the line charac- 
teristics with those of the transmitting filter into a single composite 
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network function determining the wave presented to the receiving 
bandpass filter. 

The receiving bandpass filter is necessary to exclude out-of-band noise 
and interference from the detector input. It also shapes the signal wave- 
form and can include compensation for linear in-band distortion suffered 
in transmission. Two contradictory attributes are sought in the filter — a 
narrow band to reject noise and a wide band to supply a good signal 
wave to the detector. Previous work 3 has indicated a cosine filter as a 
near optimum. 

The noise-free input to the detector will be written in the form 

V r (t) = P(t) cos (cc c t + 6) - Q(t) sin (u> c t + 0). (8) 

P(t) and Q(t) represent in-phase and quadrature signal modulation 
components respectively, which are associated with a carrier wave at 
the midband frequency o> c with specified phase 0. Such a resolution can 
always be made, even though the details in actual examples may be 
burdensome. The added noise wave at the detector input is assumed to 
be Gaussian with zero mean and can likewise be written as 

t>(0 = *(0 cos (iaj. + 0) - y(t) sin (u> c t + 0). (9) 

If v{t) represents Gaussian noise band-limited to ±2co c , x{t) and y(t) 
are also Gaussian and are band-limited to ±w„ . If the spectral density 
of »(0 is w,(u), the spectral densities of a;(0 and y{t) are given by 3 

W x (u) = Wy(hi) = W„(a> c + u) + W v (b>e — CO), | 03 I < Wc . (10) 

In general, x{t) and y{t) are dependent, with cross-spectral density 

Wxy{u) = j[Wv(0)c "CO) - W v (<a c + 0))] (11) 

and cross-correlation function expressed in terms of R v (t), the auto- 
correlation function of v(t), by 

R xy {r) = -2B.(t) sinavr. (12) 

The cross-correlation vanishes at t = 0, and hence the joint distribution 
of x(t), y(t) at any specified t is that of two independent Gaussian 
variables. 

A convenient analytical model of the detector is a multiplier with 
delayed and undelayed waves applied as inputs and with a low-pass 
filter in the output to select the difference-frequency components of the 
product. In practical systems various departures from the basic model 
may offer a more convenient realization by physical circuits. The pure 
product law can be approximated by a switched modulator, an example 
of which is shown in Fig. 6(a). Here one of the two inputs operates a 
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Fig. G — Differential detector with switched modulator. 

two-transistor reversing switch in the path of the other input. The 
second input wave is in effect multiplied by a rectangular wave with 
axis crossings determined by the first input. The same result would be 
obtained from a strictly analog multiplier if the amplitude of one input 
were severely limited beforehand as shown in Fig. 6(b). We will base our 
analysis on a true productor with no limiter, but we show in Appendix 
A that limiting one input does not affect the results in the narrow-band 
case. In this detector the noise waves associated with the delayed and 
undelayed signal inputs are correlated. The amount of correlation de- 
pends on the value of the delay. The value of the delay relative to one 
bit interval will remain as a parameter to be optimized. However, we 
will require that the delay line should be designed to have a phase shift 
equal to an odd multiple of ir/2 radians at the midband angular fre- 
quency w c . 



VI. ANALYTICAL SOLUTION 



By adding (8) and (9) we obtain the signal at the input of the detector 
E(t) = xi(t) cos co c t — yi(t) sin u c t (13) 
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where 

*i(t) = x(i) +P(0 

(14) 

vi(t) = 2/(0 + Q(t). 

The delayed signal is, 

E(t - T ) = Xi(t - t) cos (u c t - o) c t) - yi(t - r) siii (w c t - <a e r) (15) 

where 

< r < T. 

In the case in which « c t = 270°, we write for E(t - r) 

E d (t) = -xus'muct — tjidcoscoct (16) 

where we set /(J — t) A/«/ . 

The low-frequency component Et t (t) of the product E(t)E d (t) is 
given by 

2Ei f = yiXid - Xiyid = {. (17) 

The performance of our system can be studied by analyzing the prob- 
ability distribution of the quadratic form in (17). Since this is a binary 
system we only need the distribution at one point, namely "zero." 

By a relabeling of the variables in which iji , x\ , iju , and xu are 
replaced by Xi , x 2 , Xz , and .r 4 , respectively, the calculation is reduced 
to the single problem of evaluating the probability that the quadratic 
form XiXa + X&3 is negative or positive, where X\ , x 2 , x 3 and x\ are 
Gaussian random variables with equal variances a 2 and mean values 
given by : 

E{x 1 ) = Q(t), E{x 2 ) = -P(t) 

E\x 3 ] = Q(t - r), E{ Xi \ = P(t - r). 

We remark that all these average values, in general, will depend on the 

signal sequence. 

A solution of this general problem in terms of an integral has been 
found 3 when the variables are independent. The present case is more 
complicated in that with an arbitrary delay r, the noise samples become 
dependent. We must now include a nonzero covariance of .x - i and xu 
and of yi and ij\ d . We point out that a solution based on uncorrelated 
noise samples would indicate a considerably poorer performance than 
found when the actual correlation is included. This is a case in which 
noise correlation is beneficial rather than harmful. 
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It has been found possible to apply the previous solution for the 
independent case by subjecting the four dependent Gaussian variables 
to a linear transformation. The new set of variables Zi , z 2 , z a , and z* 
becomes independent while preserving the invariance of the quadratic 
form on which decisions are based, i.e., 



ZiZi + ZiZ-i = X\Xi + X&i . 

The nonzero covariances of the variables are: 

cov (xi , x-i) = cov (y,y d ) = ra 1 
cov (x 2 , xa) = —cov(x,Xd) = —r<r 2 . 



(19) 



(20) 



The value of r is the normalized autocorrelation function of the noise 
evaluated at lag time t. The autocorrelation function is the Fourier 
transform of the spectral density. In the optimum receiver design the 
filter preceding the detector has a cosine amplitude-frequency response 
reaching zero at a frequency displacement from midband equal to the 
bit rate. Hence with white Gaussian noise on the line, the spectral 
density of the noise at the detector input is a squared- or raised-cosine 
function. Such a spectral function has the property that its Fourier 
transform, the autocorrelation function, decreases to zero when the 
lag time increases from zero to T. The solution we shall give is valid 
for any value of r within the permissible range from — 1 to + 1 . 

The transformation that satisfies (19) and at the same time diago- 
nalizes the covariance matrix of the z variables is: 



Z\ = |(jci — X 2 + X-i + Xi) 

z-> = |(a;i + .t 2 + x s — Xi) 

2 3 = i(— £l + X-2 + Xi + Xi) 

Z\ = \(x\ -f X-z — Xi + Xi). 



(21) 



It can be verified that (19) is satisfied when this transformation is 
applied and also that if the correlation matrix of the x's is : 



R x = 



1 r 

1 -r 

r 1 

-r 1 



(22) 



the correlation matrix of the z's becomes: 
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/?, = 



"l + r 

1 + r 










1 - r 










1 - r 



(23) 



We now recall the previously obtained result 3 for the four uncorre- 
lated Gaussian variables z x , z% , z- A , and z 4 , with the pair 2, and 2 2 hav- 
ing equal variances <m 2 and the pair z 3 and z± having equal variances 



2 

a i : 



Prob [(ziZi + 2 2 2s ) < when SA + %& > 0] = A(p,a,b) 

= \ - -^ f exp (- P V) erf [op(l - x 2 )* - fcr]<*2 

2, zyir J-i 



where 



2 
P = 



a = 



6 = 



- 2 , - 2 
Z\ + 22 

2 (To 2 

2 t 2 4 + 2 2 2 3 
2o , i(Top 2 

2224 — 2l23 



2 



a\a»p 



In our case 



a,, 2 - (1 +r)o\ 



G\ 



= (1 - r). 



(24) 

(25) 
(26) 
(27) 

(28) 



It is shown in Appendix B that the integral in (24) can be simplified 
and reduced to a two-parameter form. We also show how the general 
asymptotic expression for large signal-to-noise ratio is derived. 

The simplified version of (24) is 

A( p ,a,6) = 1 f exp [— ^ 1 <b (29) 

H> ' 2tt Jo Ll + " cos <PJ 



where 



2 

c = 



cf = 



2a 



k°- + a 2 + 1 

[(A-' + o 8 - l) 2 + 4fc 2 ]> 
A- 2 + a 2 + 1 



(30) 
(31) 



and k = b/p. 
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The signal-to-noise ratio appears only in the parameter p. The asymp- 
totic expression for large signal-to-noise ratio is therefore obtained as 
the limiting form of (29) when c 2 becomes large. The result, as shown 
in Appendix B, is: 

AW) = Q(c4) ~ ^= exp [^j,] . (32) 

It is of particular interest in our problem to write the asymptotic ex- 
pression in the form : 

G{c,d) = - e-* 3 " 2 (33) 

P 

The values of a and /3 are : 

/b 2 -f-a 2 + 1 + [(/c 2 + a 2 - l) 2 + 4fc 2 ]' 



= 



47r»a[(fc 2 + a 2 - l) 2 -f 4fc 2 ]* 
2a 2 



(34) 



k* + a 2 + 1 + [(fc 2 + a 2 - l) 2 + 4A: 2 ]*' 
These are the general results we need. 



VII. A REPRESENTATIVE EXAMPLE 



A convenient analytical representation of a band-limited FM signal 
was first proposed by Sunde. 4 This representation approximates very 
closely the actual signals generated by practical data sets. 3 In the case 
of ideal transmission, we write for the signal at the input to the detector 



E(t) = A sin w d t sin u c t — Asi(t) cos u c t (35) 



where 



•Sl 



(0 = E (-!)•&*« - nT) (36) 



Tud = tv, T & bit interval. 

The standard pulse response g(t) must satisfy Nyquist's first criterion, 
i.e., 

g(mT) = 8 m0 . (37) 

The value of b„ is + 1 for mark and — 1 for space. 

As can be seen, the signal (35) can be synthesized by exciting a net- 
work having impulse response g(t) by a series of + or — impulses oc- 
curring at integral multiples of the bit interval T. The sinusoidal com- 
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ponents of the signal are added to the response to the impulses. When 
delay distortion is present, the impulse response g(t) is suitably modified 
by introducing a quadrature component. 
The delayed version of (35) becomes: 

EM) = A sin (u d t - u,iTb) sin ( Wf / - w c Tb) 

(38) 
- As x (t - TS) cos (uj - ta e TS) 

where b = t/T. When the noise components are added and the condition 
u e TS = 3t/2 is satisfied, the low-frequency product is 

£($) m -[A sin (co,,/ - irb) + .//(/ - TS)][Aai(t) + x(t)} 

+ [A sin uj + !j(t)][A Sl (t - TS) + x(i - TS)]. 

The function £(f) represents the signal as it appears at the output of 
the ideal low-pass filter. To obtain the information, £(t) must be sam- 
pled at integral multiples of the bit interval. The sample thus obtained 
at t = T + eT, ^ e ^ 1 is 

£(T + Te) = -{A sin [*-(! + e - 5)] + y[T(l + e - «)]} 

•{yl Sl [7'(l + 6)]+.r[r(l + e)]} 

(40) 
+ {.4 sinkd + e)] + Z/[T(1 + e)]j 

•Uai[r(l + e - 8)} + x[T(l + e - «)]}. 

As is characteristic of nonlinear detection processes, the presence of 
noise introduces dependence on signal history. The memory can be 
minimized by using a pulse spectrum which satisfies the second as well 
as the first of Nyquist's criteria, 5 i.e., one which preserves the spacing 
of transition times as well as axis crossings. To illustrate let e = 
and 5 = J. In this case (40) reduces to: 

«T) = -\A + y(T/2)\{ Sl (T) + as(T)| 

+ \!/(T)\\A Si (T/2) + x(T/2)} 



(41) 



/here 



Si 



(J) = E (-l)"fc^(* - nf) (42) 



Sl (r/2) = - }(bi - bo), SxCT 1 ) = -6, . (43) 

The memory is thus reduced to one previous symbol. It follows, there- 
fore, that 

S(T) = [x - Abti-ya - A] + y[x d - (A/2)(6, - b )\. (44) 
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The probability of error is the weighted average of the conditional 
probabilities that the sample is negative when 61 is +1 and that the 
sample is positive when 61 = —1. Since the sample depends on the 
present and immediately preceding symbols, there are four different 
cases to consider: 61 = —60 = 1; 61 = 60 = 1; &i = —60 = — 1; and 
bi — 60 — — 1. If marking and spacing symbols are equally probable, 
each case has a probability of one-fourth. We now write for the average 
probability of error 

4 
Pe = I Z Prob [(x ln x 4n + x 2n x 3n ) < 0] (45) 

n=l 

where the variables x in , x 2n , x 3n , and a- 4n are specified by Table I. 

A physical interpretation of the detection process in this case can be 
obtained by regarding the quadratic form in (44) as the scalar product 
of two vectors, e.g. : 

x&i + x 2 x 3 = u-v (46) 

where if i,j represent unit vectors along rectangular coordinate axes, 
the four possible pairs of vectors are : 



Ui = i(-A + x) + yy 

Vi = i(-A - y d ) + j(-A + x d ) 

u 2 = i(— A + x) + jy 

v 2 = i(-A - y d ) + jx d 

u 3 = i(A + 3) — yy 

v 3 = \{A + ya) + j(A + x d ) 

u 4 = \{A + x) - yy 

v 4 = i(A + y A ) + jxd . 



(47) 



(48) 



(49) 



(50) 



Occurrence of error is synonymous with a negative value for the 
scalar product of any pair of vectors, and hence is also equivalent to an 

Table I 



n 


Xl„ 


*»n 


Xln 


ft* 


1 

2 
3 

4 


x - A 
x - A 
x + A 

x + A 


y 

y 

-y 

-y 


Xd — A 

Xd 

Xd + A 

Xd 


-yd - A 

—yd - A 

yd + A 

Vd + A 
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angle greater than 90° between the vectors of a pair. We can use this 
representation to explain the beneficial effects of correlated noise 
samples. There is positive correlation between x and xa as well as 
between y and y,i . This means that Xd is likely to have the same 
sign as x and similarly for ijd and y. Inspection of the vector com- 
ponents shows that such noise values have opposite effects on the size 
of the angle between the two vectors, as shown in Fig. 7, illustrating 
(49). 

At this point we can obtain an explicit expression for the probability 
of error. From Table I, we note that the average values of the x's are 
given by Table II. From (21), the z's are given by Table III. It is seen 
that the cases of n = 3 and 4 differ from n = 1 and 2 respectively by a 
change of sign throughout. A change of sign in all the z's in (2o)-(27) 



y AND yd OF SAME SIGN HAVE OPPOSING EFFECTS IN GOING FROM 9 -»»i 


I AND Id OF SAME SIGN HAVE OPPOSING EFFECTS IN GOING FROM fl, -»8 Z 


4J x d 
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Fig. 7 — Effect of correlated noise samples on angle between vector inputs to 
detector. 
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Table 11 













n 


'in 


Xin 


■Tin 


Xin 


1 


-A 





-A 


-A 


2 


-A 








-A 


3 


A 





A 


A 


4 


A 








A 



leaves the parameters a and b unchanged. Therefore, the four probabil- 
ities of (45) reduce to two distinct ones with the multiplying factor 
changed from j to \. 

The values of the parameters for the two cases are found to be 



5A 2 



p\ 



4(1 + r)a 2 



a 2 A /l+l 

b, - - 



2 

P2 = 



02 = 



2(1 + r)a 2 



V 1 - r 



(51) 



2«rV5(l - r) 



/;•» = 0. 



As demonstrated in Ref. 3, the mean square of the signal wave on the 
line with random data is A 2 in the optimum case for white Gaussian 
noise. The noise input to the detector has a squared-cosine spectral 
density function, with a' equal to the area. Hence a" is one-half the mean- 
square value of the noise on the line in a band of width equal to twice 
the bit rate. In terms of the parameter M, defined as the ratio of average 
signal power on the line to average noise power in the bit-rate band- 
width, we have 



M = A 2 /a\ 



(52) 



Substituting the appropriate correlation value of r = £, we obtain 



Table III 



n 


Zin 


ijn 


Zan 


Z«n 


1 
2 
3 
4 


-3 A/2 
-A 
3 A/2 
A 


-A/2 


A/2 



-A/2 


A/2 



-A/2 
-A 

A/2 

A 
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for the two sets of parameters 



2 5 , M 

Pi =gM fn -y 

oi - | \/3 a 2 = V3 (53) 

o 

'- = - j/| '* = o. 

The probability of error is given by 

1\ = \ A( Pl , a, , &i) + | A( P2 , a 2 , b 2 ). (54) 

Using the computer program previously established for A(p,a,b), S. 
Habib has obtained the uppermost curve of Fig. 8. Also shown are the 
ideal performance for coherent binary detection, the ideal performance 
for noncoherent binary detection, and an experimental curve obtained 
by E. R. Day. 

It is also instructive to apply the asymptotic formulas given in (33). 
In this example, applying (53) to (34), we obtain 

3 /T . 2 

4 V Sir 5 



«i = : - A/ — , fix = 1 



(55) 



2V2 



L'tt 

Substituting these values in (30) we find 

'1 3 Q 

r» - u ' 3 1 ..~ J//3 ,-' u/3 /-r\ 

p '~W*rti e + 8V2^7 e ~8V2^7 e • (56) 

In terms of the notation introduced in Section IV, k = \. The asymptote 
is plotted in Fig. 8 and is found to be very close to the result of the 
exact calculation in the region of interest. Since the optimum coherent 
system has an error probability proportional to e~~ , the half -bit differ- 
ential system in the limit requires 10 log™ 3 or 4.8 db more signal-to- 
noise ratio than the optimum for the same performance. Of this penalty, 
3 db is accounted for by the steady-state informationless tones which 
make up half the average power of the FM signal. The remaining 1.8 db 
can be ascribed to the differential-detection scheme. 

It is also possible to decode the message by providing a full-bit delay 
at the receiver. However, the decoding has to be performed on the 
transitions. In addition, the phase shift must be such that cos u) c T = 1. 
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Fig. 8 — Error rates in binary data transmission. 

As before, if we perform the indicated multiplication, reduce to a sum 
of single-frequency terms, and reject those that contain the midband 
frequency co r , we obtain the post-detection filter output, which we now 
designate as Vi f (t). Again, to simplify the formulas we omit showing 
the functional dependence on t and use the subscript d to indicate 
values of functions at t — T. Then 



2V [f = XiXu — XiAsu — xuAsi + A SiSu 



+ A(yi — yu) sin u d t + ij\\ju — A 2 sin 2 u d t. 



(57) 



(01) 



[62) 
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Finally, we assume samples to he taken at multiples of T, which means 
that sin tad becomes sin u<mT = and 

Sl = Sl (mT) = (-Y"b m (58) 

8id - «i[(« ~ I)?'] = (-)"'^ m -i • (59) 

For the rath sample, then 

2Vi f = XiXu + IhUui - (-)'"A(b ni x ld - b m -iXi) - .4' 2 6„,6„,_i . (60) 

Let the binary information be coded in the transitions. Then, for ex- 
ample, a "1" is represented by b m = —b m -i and a "0" by b,„ = b„_i . 
The two signaling conditions are then: 

2Vi/i = A a + (-)'"b m AU i + x u ) + xix u + ywid 

= [A + (-)"7>„,.r,][.4 + (-)"'&,„»•„,] + //.//.,/ 

2V U0 = -A' 2 + ( - ) "Kit (xi - xu) + XiXu + yiyu 

= -[.1 - ( - )"'b,„.v 1 }[A + (-)'"/>,„.(•„/] + //,//„/. 

The sampling interval T is typically large enough in this case to make 
the delayed noise samples independent of the direct samples. With this 
assumption we can regard X\ , v/i , Xu , and yu as independent Gaussian 
variables. In the first condition, we set 

H -= A+ {-)'"b m .i\ (63) 

U = A + {-) m bn&u (64) 

27„ = &+ yiVid. (05) 

Then £, £</ , >/i , and ,(/i,/ are independent Gaussian variables with standard 
deviation a, where a is the rms noise voltage at the detector input, i.e., 
the rms value of either X\ or y, . The mean values of iji and yu are zero, 
and the mean values of £ and £,/ are .4 . In the second condition, we set 

g = -,1 + (-)"7wi (66) 

and obtain the same relations except that the mean value of £ becomes 
—A instead of .4. 

Correct decisions are made in the first condition if 2Vi/ is positive, 
and in the second condition if 'IV u is negative. Hence if we let 

2 = & + UWu (67) 

and designate pi(z) and pi(z) as the probability density functions of z 
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when £ has the means A and —A respectively, the probabilities of error 
in the two conditions become 



Pi = [ pi(z)dz 

v— 00 

Pi = I -pi(z)dz. 
Jo 



(68) 



(69) 



This is the same problem solved in equations (54)- (58) of a previous 
paper 3 , and the final result is found to be 

„-A*/(2«*) __ i „-MI2 ,jq\ 



P = P 1 = P 2 =I e - 



_ 1 
" 2 



= •. C 



Equation (70) shows 3 db poorer performance than ordinary differ- 
entially coherent phase modulation. The one-bit delay differential 
system with transition coding thus suffers 3.9 db penalty relative to 
ideal coherent detection in the error probability range of 10~ . 

An interesting result is obtained when the delay line at the receiver 
is allowed to become small relative to the bit interval while still main- 
taining the condition that the phase shift at the carrier frequency ta e 
equals 270°. We show in Appendix C that the performance in this case 
approaches the performance of an ideal phase differentiator. 

VIII. RESULTS WITH NO DELAY DISTORTION 

The first part of our numerical results deals with the performance of 
the differential FM detector as a function of the value of differential 
delay and the sampling instant in the absence of delay distortion on the 
line. As a preliminary, we show in Figs. 9-12 inclusive the computer 
print-outs of the eye patterns for differential-delay values of 0.2, 0.4, 
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Fig. 9 — Eye pattern for differential delay of 0.2-bit interval. 



BINARY FM DETECTION 



135 



o 
-2 



-0.5 0.5 

TIME IN BIT INTERVALS 



Fig. 10 — Eye pattern fur differential delay of 0.4-bit interval. 
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Fig. 11 — Eye pattern fur differential delay of O.G-bit interval. 
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Fig. 12 — Eye pattern for differential delay of 0.8-bit interval. 
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0.6, and 0.8 bit interval. Some of the features found have been previously 
mentioned in Sec. III. As before stated, the patterns are not to be 
interpreted as furnishing quantitative evaluation of performance. For 
the latter, we rely entirely on error rates vs signal-to-noise ratio. The 
figures illustrate what would be seen at the detector output in the 
absence of noise for the various cases and indicate the preferred sam- 
pling time for each case. 

For any pair of values of 8 and e in (40), £ is a random variable possess- 
ing a probability density determined by the additive noise and previous 
signal history. We have conserved computation time without loss of 
essential information by concentrating attention on the asymptotic 
performance when the signal-to-noise ratio is large. The corresponding 
db impairment or degradation relative to optimum binary PM as 
expressed by 10 log™ (1/x) has been computed for various values of 
8 and e and for all data sequences of length 5 bits. For a few representa- 
tive cases, we have used the exact formula over a range of signal-to-noise 
ratios to indicate the degree of approximation given by the asymptotic 
formula for typical conditions of interest. A raised-cosine pulse spectrum 
on the line has been assumed. Details of the computations are given in 
Appendix D. 

The IBM 7090 computer was programmed to evaluate the parameters 
developed in Appendix D for all 32 possible 5-bit sequences and for 
different values of 8 and e. Figs. 13-18 represent the degradation vs 
value of delay for different sampling instants across the bit. The sampling 
time is measured from the midpoint of the bit interval in the undelayed 
wave. The curves on each graph represent the various sets of sequences. 
It appears from the graphs that the 32 sequences tend to bunch into 
four distinct sets. We did not attempt to label or identify these sets, 
since the average performance is more closely determined by the se- 
quences which suffer the most. This is because the relative probability 
of error for the sets varies exponentially with the degradation. It can be 
seen that as the value of delay exceeds half a bit the performance 
degrades rapidly. 

We next examined the performance of the receiver for fixed delay line 
values and variable sampling instants. As shown in Appendix C, the 
d<p/dt detector may be regarded as a limiting form of differential delay 
as the value of delay approaches zero. Fig. 19 shows the performance 
for the worst and best sequences as the sampling instant is varied across 
the bit. It is clear from this figure that the best sampling instant is in 
the middle of the bit, i.e., midway between transitions of the detector 
output wave. It turns out that the best sequence is the sequence of all 
marks or spaces, while the worst sequence is that of reversals. 
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RATIO OF DELAY LINE TO BIT INTERVAL 

Fig. 13 — Db degradation vs differential delay with sampling at 0.005-bit 
interval. 
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Fig. 14 — Db degradation vs differential delay with sampling at 0.11-bit in- 
terval. 
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Fig. 15 — Db degradation vs differential delay with sampling at 0.21-bit 
interval. 
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Fig. 16 — Db degradation vs differential delay with sampling at 0.31-bit 
interval. 
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Fig. 17 — Db degradation vs differentia] delay with sampling at 0.41-bit 
interval. 
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Fig. 18— Db degradation vs differential delay with sampling at 0.51-bit 
interval. 
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Fig. 19 — Db degradation vs sampling time for d<p/dt detector with no delay 
distortion. 



Fig. 20 exhibits the same information as Fig. 19 with the half-bit 
delay. In this case the worst sequence is a function of the sampling 
instant. At one-fourth bit and three-fourths bit sampling points re- 
ferred to an origin at the transition points of the detector output, the 
intersymbol interference is zero, and at these instants the degradation 
is 4.8 db for all sequences. Nearer the center of the bit interval as seen 
from the detector output, the degradation is slightly greater. We never- 
theless conclude that the best sampling instant is at this center, which 
we shall refer to as "mid-bit." By sampling at this time, we obtain a 
spread of half a bit interval tolerance to sampling jitter. 

Similar curves are found when the receiver delay line is other than 
zero or a half-bit. In Fig. 21 we show the degradation of the worst 
sequence as a function of the receiver delay line. The zero and half-bit 
values are as shown on Figs. 19 and 20. The curve as shown applies to 
use of a delay line in which the phase shift is an odd multiple of 7r/2 at 
midband. We have previously shown that the signal can be recovered 
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Fig. 20 — Db degradation vs sampling time for half-bit differential-delay de- 
tector with no delay distortion. 
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Fig. 21 — Db degradation vs differential delay for most vulnerable sequence 
sampled at mid-bit. 
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with low degradation at a full-bit delay by differential encoding and use 
of a delay line with phase shift equal to a multiple of t at midband. 

The above results indicate that, neglecting sampling jitter, least 
degradation results with least receiver delay, and the dip/dt receiver is 
best. Differences are small up to a half-bit. Addition of delay distortion 
alters this conclusion. 

IX. EFFECTS OF DELAY DISTORTION 

We preface our discussion by exhibiting sets of computed eye patterns 
for cases of linear and quadratic delay distortion in the channel. The 
amount of distortion is specified by the increment in delay measured in 
bit intervals between the center and edge of the transmission band. 
Results for the dip/dt detector are shown in Fig. 22 and for the half-bit 
delay differential detector in Fig. 23. The eye is found to close for smaller 
amounts of linear delay variation than for quadratic. As before, these 
results are not to be taken as quantitative measures of performance. 
The same numbers which determine the eye traces arc also used in 
calculating error probabilities, but in a different way which precludes 
derivation of cither final result directly from the other. 

Fig. 24 shows the calculated degradation from the ideal for the dip/dt 
detector with various amounts of quadratic delay distortion measured 
in bit intervals. The results are plotted as a function of sampling time 
for the most vulnerable data sequence. There is an indication that the 
best sampling time is not at mid-bit when the delay distortion is large. 
However, the decreased tolerance to timing jitter would make such a 
shift undesirable. The effect of linear delay distortion, as exhibited in 
Fig. 25, is considerably worse. A comparison between effects of the two 
kinds of distortion is obtained in Fig. 26 by replotting the mid-bit 
sampling results of Figs. 24 and 25 as a function of delay distortion. 
As stated before, the dip/dt detector is equivalent to zero differential 
delay. 

Figs. 27 and 28 present corresponding curves for the case of a half-bit 
differential delay. The shapes are similar to those for dip/dt. There is 
slightly more tolerance to jitter, and the best sampling time is still at 
mid-bit. 

In an effort to compare various receiver delay lines, previous data arc 
cross-plotted in two ways. In Fig. 29, we in effect extend Fig. 21 to 
show how performance is affected by amount of differential delay when 
various amounts of specified linear and quadratic delay distortion are 
present in the channel. The curve for zero delay distortion replotted 
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Fig. 24 — Db degradation of d<p/dt detector vs sampling time for quadratic 
delay distortion. 

from Fig. 21 shows the dip/dt receiver as best for a constant-delay channel. 
For quadratic delay distortion of two bit intervals or linear delay distor- 
tion of one bit interval, there is slightly less degradation at a half-bit of 
differential delay. For larger amounts of delay distortion the advantage 
from more receiver delay is greater. 

The same conclusions can be drawn from Fig. 30, in which the abscissas 
are amounts of quadratic delay distortion and the curves are drawn for 
specified values of differential delay. The curves cross over from the 
condition of a preference for least differential delay with low delay 
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SAMPLING TIME IN FRACTION OF BIT 



Fig. 25 — Db degradation of d<p/dt detector vs sampling time for linear delay 
distortion. 



distortion to a preference for greatest delay with high delay distortion. 
The exact choice of delay line depends on several factors. Say the maxi- 
mum delay distortion is not greater than three bit intervals, then the 
half-bit delay is best over most of the range and only slightly worse 
over the rest. If the delay distortion is never more than one bit interval, 
a dtp/dt detector is best, while if very high delay distortions are encoun- 
tered, and fairly high degradation is permissible at lower values, a delay 
line of 0.7 bit interval would be better. It was also found that the longer 
delay lines provide more tolerance to jitter in the sampling time. 
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Fig. 28 — Db degradation with half-bit differential delay vs linear and quad- 
ratic delay distortion. 

The asymptotic values of db impairment at large signal-to-noise 
ratio have been used throughout our discussion of Figs. 24-30 as a 
measure of performance for conditions of practical interest. As a check 
on the validity of this concept, complete curves of error probability vs 
signal-to-noise ratio have been computed from the exact formulas in 
representative cases. These curves are shown in Fig. 31 together with 
the asymptotic approximations. It appears that the latter are sufficient 
for most engineering applications. 



APPENDIX A 



Assume that a limiter is inserted in the undelayed input to the mul- 
tiplier as shown in Fig. 6(b). The input to the limiter is then given by 
(13), which can also be written in the equivalent form 



E(t) = R(t) cos [ujl + ip(t)] 

R(t) = [*!*(<) + 2/i 2 (0]» ^ 
tan <p(t) = yi{t)/x x {t). 



(71) 
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Fig. 29 — Db degradation vs differential delay for various amounts of linear 
and quadratic delay distortion. Curves are designated "L" for linear and "Q" 
for (juadratic, followed by number of bit intervals of delay distortion. 
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Fig. 30 — Db degradation vs quadratic delay distortion for various amounts 
of differential delay. 
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If the limiter is ideal, its output E L (t) is a positive constant E when 
the value of the cosine function is positive and is a negative constant 
— En when the cosine is negative. That is, if z = u c t + <p(t), 

\ Eq, cos z > 

EM = 

[-E , cos z < , 72 , 

_ 4Eoy (-)" cos (2w + l)z 
ir „=o 2n + 1 

Multiplying E L (t) by the delayed signal and substituting cos co c t = 0, 
sin u c t = —1, gives: 

E L (t)E(t - t) = 2^ ^ , , 

7T „=o 2n + 1 

■ [sin \2nu,.t 4- (2n 4- 1 VU) - *>0 - r)] (73) 

- sin [2(n + 1 W + (2n + 1)*>(0 

+?(/ - r)]j. 

If it is possible by low-pass filtering to accept the term 

sin [<p(t) - <p(t - t)} 

while rejecting the next higher-frequency terms sin [2a) c l + 3<p(/) — <p 
(t — t)] and sin [2a>,/ + <p(t) + <p(t — t)], the filtered response becomes: 

Vt At) = 2I ' J ° R(i - T) 8 in bit) - pit - t)] 



IT 

2E t) R(l - t 



(sin <p cos <p,i — cos ip sin tpi) (74) 

7T 

_ 2EpR(t — r)(yixi d — Xiy u ) 
*{x? + y?)K*rf + 2/w 2 )* 

The only term in (74) which can have both plus and minus signs is 
yiXid — xijjid , which is the same term found to be the basis of binary 
decisions in (17) for the case in which a pure product was taken. The 
switched modulator therefore gives the same error performance as the 
productor if there is sufficient frequency separation between the desired 
low-frequency output and the sidebands on 2a> e . 

An ideal limiter was assumed for simplicity in the argument just 
given, but the equivalence can be proved for a wide class of nonlinear 
devices in one of the two paths. It is sufficient that the output of the 
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device is of the same polarity as the input and can be expanded as a 
Fourier series in terms of the input values. We can then write 

E L {t) = F(R cos z) 

= tt + 1~> * n cos nz 

Z n— 1 

F„ = - I F(R cos z) cos nz dz. (70) 

T J-ir/2 

Note that 7? varies with time but cannot be negative. 
In particular 

I ,.x/2 , ^.31/2 

F, = - / F(R cos z) cos z dz + - / F(fl cos z) cos 2 dz. (77) 

T J- rll T J- til 

In the first integral cos z is positive and hence F(R cos 2) is positive. 
In the second integral cos 2 is negative, F(R cos 2) is negative, and the 
product is positive. It follows that F\ is positive. 

If we repeat the calculation leading to (74) with the Fourier series 
(75) replacing (72), we find that 2E /ir is replaced by Fi , and since 
F\ is positive, the conclusion is the same. It will however be more diffi- 
cult to isolate the desired low-frequency term when Fn and F 2 are present, 
since these coefficients lead to components centered about w r . It is 
therefore preferable that the nonlinear function F have odd symmetry 
about the origin, in which case the even-order coefficients vanish. 

The more general argument is useful in showing that the switched 
modulator does not have to be perfect. Note that the equivalence is 
destroyed if limiters are inserted in both paths. There would then be 
sidebands on harmonics in both inputs to the modulator, and the beats 
between these sidebands would generate additional components in the 
low-frequency band. 

APPENDIX B 

Simplification of the Error Probability Formula, and Determination of Its 
Asymptotic Form 

First, replace the parameter b by 

/, = b - = 1^-^. (78) 
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The expression to be evaluated is now 



Ai(p,a,fc) = g - ^T7^ I e- p2xi erf [pUvT^T 2 - kx)]dx. 
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(79) 



Substitute px = x' and replace the error function by its definition as an 
integral. The result is 

1 ! ' ' e- {xHv) dxdy. (80) 



Ai(p,a,A-) / / 

2 7T J-p JO 



A graph of the region of integration in the xy-p\ane is shown in Fig. 32. 
From this figure, we deduce that by a change to polar coordinates, we 
obtain the equivalent expression 




Fig. 32 — Region <>f Integratiun for evaluation of double integral in (80). 
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.x— arc tank /-AC) 



. , «x— arc tans fJlKPI 

Ai(p,a,k) - - — - J dd e" r V dr 

L IT Jo J 

1 *V r—plcUB S 

- - de e~ r \dr (81) 

IT J r— arc tan k J 

1 »— arc tan k /-p/cos $ 

— - I dd I e~ r r dr 

T Jo Jfi(B) 

where 

fi(8) = Ap[(sm d + k cos ef + a cos 2 d]~ m . (82) 

The integration with respect to r can be performed, and after some sub- 
sequent combining of terms, we find 

A,(p,a,fc) = ^ f ' exp ["- . . — - t ^-n 2 nil *• (83 > 

2t Jo L (sin + k cos 0) 2 4- a 2 cos 2 0J 

By the substitution 20 = <p, (83) transforms to 

1 r 2 * 

Ai( p,a,k) = — / 

47T • u 

' eXP L~/o 2 + a 2 4- 1 4- (fc 2 + a 2 - 1) cos <p - 2k sin <J ^' 
We then note that 
(A: 2 -\- a — 1 ) cos (p — 2k sin <p 

= [(k 2 + a 2 - l) 2 + 4A: 2 ]* cos L + arc tan fe2 + 2 q t _ X 



(84) 



(85) 



Taking advantage of the fact that the range of integration is one com- 
plete period of the integrand in <p, we can replace the sum of <p and a 
constant angle by a new variable without changing the limits. Noting 
furthermore that the integrand then becomes an even function of the 
variable of integration, we obtain finally 

Ai (p,a,k) = G(c,d) 

1 r f c 2 1. (86) 

= 2TJo eXP L" l+^cosJ ^ 

where c and a have the values given by (30) and (31). 

Applying the method of steepest descents to the case in which c is 
large, we write 
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G(c,d) = i- C e- c2 * w d<p (87) 

lir Jo 

$(#) = (i + d 2 cos v)- 1 

&U) = d\l + d 2 cos <p)~ 2 sin ^c 
v ^ s (88) 

3>'(<po) = at <po = or 7r 

$"(«5o) = d 2 (l + rf 2 cos v?o) -2 cos (p > for vr (1 = 0. 
Hence, we set <p = and approximate $(<?) by 

<p(<p) = <S>(«po) + 5*"(vo)(^ — Vo) + • • • 

= (i + rf 2 ) - I + ^ 2 (i + rf 2 rV+ •••• 

If c 2 » 1, 



(89) 



1 -f- r/ 2 



2CC/V27T 



= exp 



(-Thy 



(90) 



APPENDIX C 



Limit of Performance as Delay Approaches Zero 

As the delay t is made small, the delayed variables in (17) can be 
expressed in terms of the undelayed ones by the linear approximations, 

x id tt Xl — TXi 

(91) 
V\d ^ 2/i - rill , 

where the dot signifies the derivative with respect to time. In the limit 
we then find 

2E lf K r{x x y x - y^) . (92) 

The amplitude of the detected signal approaches zero as the delay is 
made small, but if t is not actually zero, the lack of output can be com- 
pensated by linear amplification. Hence in the absence of imperfections 
other than additive noise in the channel, binary decisions are made on 
the basis of the sign of xyi — y&i . 

In the case of a dtp/dt detector, binary decisions are made on the 
basis of the sign of 
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^-^togt- ^-f* . (93) 

dt dt Xi Xi 2 + ?/i 2 

Since X\ + ?/i 2 cannot change sign, the decisions are actually made on 
the sign of xiy\ — yi±i . We conclude that the binary error rates of the 
two systems must approach equality as t approaches zero. 

APPENDIX D 

Computational Details 

For the pulse response (36) we use the impulse response corresponding 
to a raised-cosine spectrum, namely 

sin (y') 

9(0 = /oA r T^Tvn- (94) 



•©[■-©•]■ 



The sampled value of the pulse train at t = T + eT is then given by 

«i[T(l + 01 

V /■ i^i. sin2x[l - 8 2 - n] _ o, s v ( 95 > 

= V (_1) 6 " Ml - & - n)[l - 4(1 - 5 2 - n)1 = * M 

and the delayed version after r seconds by 

Sl [T(l + «-!)] 

v , n », sin 2tt[1 -ji-w] «/.x (96) 

» 2rr(l — 5i — tt)[l — 4(1 — 81 — ?i) 2 \ 

where 

8 = t/T, 8 - e = 81, Si- S = Si. (97) 

In terms of the two new variables Si , 5 2 , the noise correlation is 
given by 

/.\ sin 2tt(5i - 82) , . 

• W = 2t(5< _ m _ 4(5i _ fc)1) • (98) 

The observed voltage or current £ can now be written as a function of 
81 and 82 in the following form : 

£(5i , 5 2 ; = [Asi°(8 2 ) 4- x][-A sin ttSi - y d ] 

(99) 
+ [A sin tt5 2 + y][A*i M + «d • 



p = 



&A* 1 + rfa - 82) 

ZiZj + Z2Z3 

!i 2 + U 



,/l 


+ 


rfa - 


- 5 2 ) 


\ 1 


— 


r(«, ■ 


- h) 


/l 


+ 


rfa - 


- h) 



(100) 
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In the calculation of eye patterns, the noise samples x, x* , y, and ya 
are omitted. By dropping these terms in (99) and substituting suc- 
cessively e = and e = 8, we verify that the samples at these instants 
depend on only one value in the data sequence. This confirms the state- 
ments made in Section III in the discussion of eye patterns. When the 
noise samples are inserted, the entire previous signal history exerts its 
effect at all sampling instants, and in fact the preferred sampling instant 
from the standpoint of low error rate and tolerance to jitter in sampling 
time is not necessarily the one in which intersymbol interference vanishes 
in the absence of noise. 

The probability of error for a particular sequence is : 

P c = § Pr {[AsSfa,) + x][-A sin tt5i - y,\ 

+ [.4 sin tt5 2 + y][Asi + fa) + M < 0| 
+ \ Pr |[^sr(5 2 ) + x][-A sin x5 x - y d ] 
+ [A sin tt5.. + y][Asffa) + x d ] ^ 0} 
where 

_ Vr _ „, sin 2t[1 - n - n] < 101 

fo "2tt(1 -m-m)[1 -4(1 - M - n) 2 ] ± r(/x) " 

Applying the transformation in (21) we obtain: 

z* = (A/2)[si ± (5 2 ) - sin irS 2 + 8* fa) - sin ^ 
l± = (A/2)[s 1 ± (5 2 ) + sin tt5 2 + s^fa) + sin ir5i] 
it = (A/2)[-stfa) + sin t8-> + at fa) - sin t5j] 
zt = (A/2)[stfa) + sin tt5 2 - 8* fa) - sin ttSi] . 

The required parameters in (25)-(27) are 
M I, 2 + zV 
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(103) 

(104) 



^^v rxTv (io5) 

2l 2 + Z2 2 K 1 — r(5i — 5 2 ) 
As pointed out before, the asymptotic degradation is given by 

degradation = 10 log,„ (l/«)db (106) 
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where 

fc 2 + q 2 + l + [(fc 2 + a 2 - l) 2 + 4fc 2 ]* ( . 

K = W ~~ (10?) 

and 

P = M Pl . (108) 

When delay distortion is present in the channel, a convolution is 
performed to evaluate the resulting in-phase and quadrature components 
of the noise-free input to the detector. It is convenient to combine the 
assumed delay distortion with the equivalent amplitude characteristic 
arising from signal pulse shaping, sending filter, transmission line, and 
receiving filter to form a single complex transmittance function 

ff(u) = A(<a)e Mu) . (109) 

The impulse response of the medium is then 

hit) = I H(w)e ial d«/2ir 

J- oo 

= f A(o>)e ma)+al1 du>/2ir (110) 

= / A(w) cos [d(co) + ut] doi/ir. 
o 

Let A(u) = B(co — u) c ) = B(v), 0(co) = <p(u} — Uc) + v?o = <p(v) + vo , 
and decompose (110) into the in-phase and quadrature components 



where 



h(t) = I B(v) cos [<p(v) + vt + co c t + <po] dv/ir 

■'-wo 

= hi(t) cos u) c t — h 2 (t) sin u c t, 
hi(t) = I B(v) cos [<p(v) + vt + <po\ dv/ir 

J-ttQ 

hi(t) = I B(v) sin [<p(v) + vt + po] dv/ir. 

J— tin 



(111) 



(112) 



Since the medium is assumed to be linear, the signal input to the 
detector is the convolution of the input and the impulse response, 

V 1 (t) = [ V(t- T )h(r)dr, (113) 

J— 00 
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where 7(0 can be written in the form: 

7(0 = Po(0 cos aj - Q (0 sin "J- (H4) 

Then from (111), 

7i(0 - ( [P«(t - t) COS («J - o> c r) 

- Qo(< - t) sin (ujt - w e T>][Ai(r) cos co c t - fci(r) sin co c r]f/T. 

Dropping the double-frequency terms, we obtain for the input to the 
detector 

7,(0 = P(t) cos tad - Q(t) Bin taj (116) 



where 



P(t) = ^ [P (t - r)Ai(r) - Qo(t - T)h 2 (r)}dT/2 

J— oo 

Q(t) = [ [P (t - r)hM + QoU - r)hi(T)\dr/2. 

J— 00 



(117) 



REFERENCES 

1. Davey, J. R., Signal Space Diagrams, B.S.T.J., 43, Nov., 1964, pp. 2973-2983. 

2. Brand, S., and Carter, C. W., A 1,650-Bit-Per-Second Data System for Use over 

the Switched Telephone Network, A.I.E.E. Trans., Pt. I, Coram, and Elect., 
80, Jan., 1902, pp. 652-661. 
3 Bennett, W. R., and Salz, J., Binary Data Transmission by FM Over a Real 
Channel, B.S.T.J., 42, Sept., 1963", pp. 2387-2420. 

4. Snnde, E. 1)., Ideal Pulses Transmitted by AM and FM, B.S.T.J., 38, Nov., 

1959, pp. 1357-1420. 1TI1T1 

5. Nyquist, H., Certain Topics in Telegraph Transmission, Trans., A.I.E.E 

47, pp. 617-644, Apr., 1928. 



